Quantum waveguide theory of Andreev spectroscopy in multiband superconductors: 

the case of Fe-pnictides 
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The problem of Andreev reflection between a normal metal and a multiband superconductor is ad- 
dressed. The appropriate matching conditions for the wave function at the interface are established 
on the basis of an extension of quantum waveguide theory to these systems. Interference effects 
between different bands of the superconductor manifest themselves in the conductance and the case 
of FeAs superconductors is specifically considered, in the framework of a recently proposed effective 
two-band model, in the sign-reversed s-wave pairing scenario. Resonant transmission through sur- 
face Andreev bound states is found as well as destructive interference effects that produce zeros in 
the conductance at normal incidence. Both these effects occur at nonzero bias voltage. 
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I. INTRODUCTION 

Electronic scattering at the interface between a nor- 
mal metal (N) and a superconductor has been used 
as a probe to investigate the electronic properties of 
superconductors 1 "^ and, more recently, FeAs supercon- 
ductors(FAS), leading, in the latter case, to different con- 
clusions regarding the pairing symmetry 3,4 . As compared 
to conventional and high-T c materials, the recently dis- 
covered FeAs based superconductors have a more com- 
plex band structure, with a Fermi surface (FS) consist- 
ing of four sheets, two of them hole-like, and the other 
two electron-like^^. S-wave, d-wave and p-wave pair- 
ing scenarios have been proposed to describe the super- 
conducting state d 10 ' 11 ' 12 . One of the suggested pairing 
scenarios is the so-called sign-reversed s-wave state (s ± - 
state), where the gap function has opposite signs in the 
hole-like and the electron- like sheets of the FS. Since this 
is a novel possibility, it deserves some theoretical develop- 
ment. A recent experiment seems to confirm this pairing 
scenario in a 122 compound 13 . 

Blonder et afl^ devised a theory for Andreev scat- 
tering in isotropic s-wave superconductors which has 
been later generalized to unconventional (anisotropic) 
superconductors^. These theories apply to one band su- 
perconductors. In the case of multiband superconductors 
(MBS), such as FAS and heavy- fermion compounds, the 
bands are usually treated as separate conduction chan- 
nels with (classically) additive conductances^, like par- 
alel resistors, thereby neglecting the quantum mechanical 
nature of the scattering problem at the interface, where 
interference effects between the transmitted waves in dif- 
ferent bands of the MBS are expected. Such interference 
effects will lead to new features in the conductance. 

We are thus posed the problem of finding the wave 
function for the scattering state of an incident particle 
from a one-band metal which is transmitted through two 
or more bands inside the superconductor. The splitting 
of the incident electron's probability amplitude among 
several conduction channels is the same quantum me- 



chanical problem as in a quantum waveguide. Thus, in 
order to derive the appropriate matching conditions for 
the wave function at the interface, we need to make an 
extension of quantum waveguide theory. 

Applying to the case of FAS, we obtain the differential 
conductance curves vs bias voltage and explicitly show 
the emergence of Andreev bound states (ABS) in the s"* 1 - 
state scenario, as a manifestation of interference effects 
between the bands, unlike the usual ABS in one-band 
superconductors. An unusual feature of the ABS is that 
they occur at a finite energy above the Fermi level and 
disperse with the electron's transverse momentum. On 
the other hand, interference effects may also suppress the 
conductance at certain energies. 



II. QUANTUM WAVEGUIDE THEORY 

The splitting of the incident electron's probability am- 
plitude among several conduction channels is the same 
quantum mechanical problem as in a quantum waveg- 
uide. In a quantum waveguide, three one-dimensional 
conductors intercept at one point (see Figure [J) 17 . The 
wavefunction for a particle must be continuous and 
single- valued at the circuit node O, implying that 



fl>( Xl -> O) = ^(x 2 -> O) = V(z 3 -» O) , 



(1) 



where x\, X2,Xs are coordinates along branches 1, 2 and 
3 respectively. The (probability and charge) current con- 
servation at the node is guaranteed by the "quantum 
Kirchhoff" lawil 



E 



1 dipjxj O) 



= 0, 



(2) 



where the coordinates Xj (J = 1, 2, 3) must be all of them 
directed to (or away from) the node O and rrij denotes 
the particle's effective mass in branch j. 

A simple one-dimensional version of the N/MBS inter- 
face is a tight-binding chain which has a bifurcation at 
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some point, as shown in Figure [2] We further assume 
the sites in branch 1 to be coupled to branch 2 through a 
hybridization operator, V. An integer n labels the two- 
atom unit cell along the chain. 



(3) 
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(2) 
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FIG. 1: Three branches of a waveguide with a node at O. 
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FIG. 2: Tight-binding waveguide with three branches. In 
branches 1 and 2 there is electron hopping along (t,t r ) and 
perpendicular (V) to the chains. An integer n labels unit 
cells along the chain. 

Let \n,j) denote the site in cell n of chain j. Then, 
the incoming particle in branch 3 with wavevector p is 
described by the wavefunction ipi nc (n) = e lpn + b e~ lpn , 
where b denotes the reflection amplitude. 

If chains 1 and 2 were decoupled, a Bloch state in chain 
j would have momentum k and energy tj{k). But now 
suppose that an operator V hybridizes Bloch states in 
the two chains. The Hamiltonian matrix for the coupled 
chains 1+2, -ffi+2, has an off-diagonal element, V(k), and 
its eigenstates follow from the eigenproblem: 



ei(fc) V(k) 
V(k) e 2 (fc) 



= £ 



(3) 



which yields two bands, £±(k), so that a Bloch state in 
the coupled chains 1+2 has the form 



4>{n) 



Oik 

Pk 



The eigenvector components, a and f3, denote the wave- 
function projections on branches 1 and 2, respectively. 
The wavefunction for the transmitted particle in chains 
1+2 reads, for n > 0, 



condition (JTJ) and by considering that the node is reached 
by formally taking n — ► 0: 



1 + b = Ca k + Da k , = Cp k + Dfa 



(5) 



In order to write Kirchhoff rule, we use the following 
expression for the probability current: 

j(r) = Re {^(r)(dH/dk)ip{r)\ (6) 

where the Hamiltonian is written in momentum space 
and the operator k = — iV in the continuum limit. In 
the tight-binding problem above, V reduces to d/dn, and 
the Hamiltonian H is just the scalar dispersion e(p) in 
branch 3, or the Hamiltonian matrix Hi + 2 m equation 
Q in branches 1+2. If we write the Kirchhoff rule as 
the following relation between the wavefunctions at the 
circuit node: 



de 

T7T yinc 
Op 



(1,1) 



8Hi 



+2 



dk 



(7) 



then, it can easily be checked that the current j(n) is 
conserved at the node, by virtue of equation ([5]). The 
left multiplication by (1,1) gives the sum of the currents 
through branches 1 and 2. Equation ([7]) reads: 



-^(1-b) 

m n 



C(l,l) 
£(1,1) 



8H 



1+2 



dk 

dHi+2 
dk 1 



Oik 
Oik' 

(3k' 



(8) 



where the effective mass m n is defined as the ratio 
between the momentum, p, and the group velocity, 
de(p)/dp. The three equations ([5]) and ((SJ uniquely de- 
termine the amplitudes 6, C, D. 

The generalization to two spatial dimensions is 
straightforward: the chain in figure [5] may be identified 
with the x direction and is repeated identically in the 
perpendicular (y) direction. The unit cell label and mo- 
mentum become two-dimensional, n and fc, respectively. 
The interface is attained as n x — > 0, the transverse mo- 
mentum component, k y , is conserved. In equation ((SJ) k 
(k') is replaced by k (fe') and the Kirchhoff rule (jHJ) is 
replaced with: 



^-(1-6) 

m„ 



C(l,l) 



£(1,1) 







dk x 1 


v Pk 


dHi+2 
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dk' x 





(9) 



ensuring the conservation of the longitudinal current j x 
at the node. 



C 



Oik 

f3k 



D 



Oi k ' 

Pk' 



ik'r 



(4) 



where the momenta satisfy the energy conservation con- 
dition £-(k) = £+{k') = e(p). We now join the wave- 
function in branch 3 with that in branches 1+2 applying 



III. MODEL FOR A FE-PNICTIDE 
SUPERCONDUCTOR 

A recent tight-binding models for the FAS band struc- 
ture assumes two orbitals per unit cell, d xz and d yz . The 
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Hamiltonian matrix is 



H(k) 



'■iy 



where \i denotes the chemical potential and 



(10) 



-2ti cos(k x ) — 2t 2 cos(ky) — 4t3 cos(k x ) cos(k y ) 
-2t 2 cos(k x ) — 2t\ cos(ky) — At% cos(k x ) cos(k y ) 
-4t4 s'm(k x ) s'm(ky) . (11) 



This is analogous to branches 1 and 2 of the waveguide 
above, with the matrix element e xy now playing the role 
of the hybridization V(k) between the branches 1 and 2 
and s x ( y )(k) playing the role of £i( 2 j(k). The parameter 
choice ti = —1, t 2 = 1.3 t 3 = t 4 = —0.85, [i = 1.45 repro- 
duces the FAS band structure^. In the unfolded Bril- 
louin Zone (BZ), the Fermi surface obtained from (flQ|) 
has two electron pockets, centered at (0, ±7r) and (±7r, 0), 
and two hole pockets, centered at (0,0) and (ir, n). 

We assume the edge of the superconductor lying along 
the y direction. Then, an incident electron on the in- 
terface with small p y is transmitted through two Fermi 
surface pockets: the electron pocket ("e FS") and the 
hole pocket ( "h FS"). See Figure [H We here work 
out the Andreev reflection problem in a FS consisting 
of just one hole and one electron pocket. The generaliza- 
tion of the theory to the four pocket FS in the reduced 
BZ or to a model with more atoms per unit ceh£ ' 19 i 20 
is straightforward. We shall concentrate below on the 
s ± -state scenario for superconductivity that has recently 
been suggested^, and show that it produces ABS as a 
consequence of interference between transmitted waves 
in the two FS pockets. 

An elementary excitation in the bulk superconductor 
with wavevector k has the wavcfunction: 



/ «k«k 



0k (r) 



ikr 




(12) 



where the coherence factors Uk,fk, denoting the ampli- 
tudes of the particle and hole components, respectively, 
obey the Bogolubov-deGennes equations^: 



H(k) 
A 



A 



I «kt»k 

UkA 

V «k/?k 



= £ 



I «kQk 

"kAc 

V w k /?k 



(13) 



with A = A(k)diag(l, 1). The superconducting gap 
A(k) is assumed to take on different values, A/j(k) and 
A e (q), in the h and e FS, respectively. In the s ± -state 
scenario, A e (k) and A^(k) have opposite signgiSI. 

The quasi-particle has a transverse momentum Hp y 
which is conserved. The incident particle from the nor- 
mal metal has momentum p + = h(p + ,p y ) and the An- 
dreev reflected hole has momentum p = h(p~,p y ). 



The transmitted particle (hole) in the superconductor's e 
band has momentum q + = h(q + ,p y ) [q~ = h(— q~ ,p y )]; 
but the transmitted particle (hole) in the superconduc- 
tor's h band has momentum = K(—k~,p y ) [k + = 
h(k + ,p y )] because the effective mass, m/j, of the h FS is 
negative and transmitted particles/holes must have pos- 
itive group velocity. See Figure [3J The wavefunction for 




FIG. 3: Schematic representation of the Fermi surfaces of the 
normal metal (left), and the superconductor's h band (mid- 
dle) and e band (right). 

a scattering state with transverse momentum hp y can be 
written as: 

tf (r) = e ip y y \^ N {x)d{-x) + ij>s{x)9{x)] 

where 9{x) denotes the Heaviside function. The wave 
function in the normal single-band metal has both parti- 
cle (u) and hole (v) components: 



4> N {x < 0) = 



a ip X 

(14) 

where a is the Andreev reflection amplitude. Near the 

Fermi level, p + ~ p~ ~ Pf\J^ — (Pv/pf) 2 , where hpp 
denotes the Fermi momentum in the normal metal, which 
has Fermi velocity vf — fipp/m n . The transmitted 
quasi-particle into the superconductor is a linear super- 
position of Bloch states of the form (fl"2"|) m the two bands: 



e lp y v ips(x > 0) 



C0 k +(r)+^k- (r) 
£<V(r) + F<Mr) ( 15 ) 



We now apply the waveguide matching conditions, at 
x = 0, between ([H)) and (fTB")) . to the u and v components 
of the wave function, separately. The condition for the 
wave function to be single valued at the node reads as: 

1 + b = Cu k +a k + + Duk-cik- + Euq+aq+ + Fu^-a^- , 

1 + b = Cu k +/?k+ + Dll k -f3 k - + i?U q +/3 q + + FlLq-flq- , 

a = Ct>k+cYk+ + -D^k-evk- + -Ev q +a! q + + Fv^-a^- , 
a = Cv k +f3 k + + Dv k -f3 k - + Ev q +(5 q + + -Fw q -/3 q - . 



(16) 



By solving the system (|16l) by the determinant method, 
the amplitudes C, D, E, F can be expressed as functions 
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of a and 6, as: 



c = 


(1- 


Vb)Ti - 


har 2 




A 




D = 


(1- 


V b)T 3 - 


- aT 4 




A 




E = 


(1- 


f &)r 5 - 


h ar 6 




A 




F = 


(1- 


f- b)T 7 - 


h ar 8 




A 





(17) 
(18) 
(19) 
(20) 

where A is the determinant of the system (|16|) and reads: 



A 



u'_a'_ 



u+(3+ u_l3_ u' + (3' + u'_{S'_ 
u+a+ w_a_ ?/a/i v'_a'_ 
V+/3+ t/^ O- 



(21) 



and the coefficients I\ are obtained from Cramer's rule. 
We now define: 



e(k) = ^(i,i) 



$(k) 



(1,1) 




dkx dk x 
for this case as: 



P- 

and write condition §^ 

1-6 = C9(k+) + L>9(k_) + £6(q+) + F0(q_ 
a = C$(k + ) + L>$(k_) + J B$(q+)+ J F 1 $(q„ 



(24) 



In order to simulate interface disorder, a potential bar- 
rier U5(x — e) is assumed in the normal metal (e < 0) and 
the limit e — > 0~ is taken22. We now show that the effect 
of the barrier amounts to making the replacement: 

1 - 6 -> 1 - 6 - 2iZ{\ + b)p F /p + (25) 
a -> a(l - 2iZp F /p + ) (26) 

on the right-hand side of equation (|24|) . where the dimen- 
sionless barrier parameter— Z = U/Hvf. To see this, we 
write the wave function in the normal single-band metal 
with both particle and hole components: 



iPn(x < e) 



e lp+x +b 



and 



(27) 



<'+< ■>■■■ Hi = n [ l Ae tp+x +p( I \, 



+ 7 







(28) 



The matching conditions at x = e < give the equations: 

TjjN(e-) = iP N (e + ), (29) 



2w„ 

Taking the limit e 

a + /3 
a — (3 

7 + <5 
7 — <5 



[V4(e + )-Vw(0] +^iv(e) =0. (30) 
0~ we obtain: 

(31) 
(32) 

= a, (33) 

(34) 



1 + 6, 

2m n U 
ih 2 p+ 



(1 + 6) + 1-6, 



2m n U 



ih 2 p 



l a. 



The waveguide matching conditions must be applied be- 
tween (|28l) and (fl"5|) at x = 0. But equations ([31]) through 
(fM)) allow the elimination of the amplitudes a, /3, 7, 8 al- 
together, finally showing that the replacements (|25|) - (f26|) 
have to be done in equation (|2~4"]) . 

The values of the Andreev and normal reflection am- 
plitudes, a and 6, can be obtained by solving the linear 
system (fl"6)) and (l24l) . Introducing 



C11 

C12 
C21 

C22 



r 1 e(k+) + r 3 e(k_) 
r 2 e(k+) + r 4 e(k_) 
r!$(k + ) + r 3 $(k_) 
r 2 $(k+) + r 4 $(k_) 



-r 5 e( q+ ) 
-r 6 e( q+ ) 
r 5 $(q+)- 
■r 6 $( q+ )- 



we obtain: 



2C21M 



and 



l + 2iZ 



\-2iZ- 



\-2iZ 



-r 7 e( q _) 
-r 8 e( q _) 

T 7 $(q_) 

r 8 $( q _), 

(35) 



(36) 



C12C2 



1 o,' 7 C22 1 C12C2 



1 + 2iZ 



1 



%• 7 C22 \ 1 C12C2 



(37) 



The contribution of this scattering state to the differen- 
tial conductance is given by: 



= 1 



\bf 



(38) 



The normal state conductance, g n = 1 — |6„| 2 , is ob- 
tained when A(k) = 0. Experimentally, the integral of g s 
(or g n ) over the transverse momenta of the incident elec- 
trons, as — J g s dp y (or cr^ = J g n dp y ) is measured 15 . 
We define the integrated relative differential conductance 
as a s /a N . 



IV. RESULTS 

The conductances as and g s {£) have been calculated 
from the above theory using the model (|10p - (fTTj) . Wc 
discuss the results below. In the calculations, the normal 
metal is assumed to have Fermi wavevector pp = 7r and 
velocity vf — 1-83. 
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When A = the reflection amplitudes, hence g s , be- 
come independent of the barrier parameter Z and this 
is precisely the condition for the occurrence of Andreev 
bound state a 15 ' 23 i 24 . Figure 2] shows the energy of the 
Andreev bound state as function of the transverse mo- 
mentum. Contrary to the usual case of single band 
non-conventional superconductors, the ABS energy is 
nonzero. It has a non monotonic dependence on p y and, 
for p y = 0, it coincides with min(|A^|, |A e |). The disper- 
sion of the ABS energy is in qualitative agreement with 
the results of Re&.. 

A = -A = 0.02 




FIG. 4: Energy of the surface Andreev bound state as function 
of the transverse momentum. A e = — A^ = 0.02. 

Figure [5] shows the conductance g s as function of in- 
cident electron energy above the Fermi level, for a clean 
[Z = 0) interface in the case where |A^| < |A e |. When 
the transverse momentum increases, g s becomes more 
strongly peaked near the energy of the ABS. In the case 
where |A/j| > |A e | there is a destructive interference ef- 
fect leading to a zero, at normal incidence, in the con- 
ductance, as shown in Figure [Sj 
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FIG. 6: Conductance g s as function of incident electron en- 
ergy for three different values of the transverse momentum, 
for a clean (Z = 0) interface. A e = 0.01, A h = -0.025. 

The effect of interface disorder is shown in Figures [3 
|S] and [HI Under increasing disorder, conductance peaks 
appear closer to the energy of ABS. The conductance g s 
is independent of the disorder parameter, Z, at the ABS 
energy, therefore, all the conductance curves g s (E), for 
different Z values, intercept at the same point, as shown 
in Figure [7] As disorder increases the conductance tends 
to decrease and, therefore, the normal state conductance 
decreases as Z increases. Since at the ABS the conduc- 
tance g s is independent of the disorder parameter, Z, 
peaks appear in the relative conductance at the ABS, 
which become more pronounced as Z gets larger. The 
relative conductance g s /g n is plotted for two ratios of the 
gap parameters in Figures [5] and O illustrating the posi- 
tions of the peaks. The case when A^ = — A e is shown 
in Figure [HI It is also seen, once again, that destructive 
interference effects in the superconductor between the e 
and h bands cause the conductance to vanish at normal 
incidence. 
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FIG. 5: Conductance g 3 as function of incident electron en- 
ergy for three different values of the transverse momentum, 
for a clean (Z = 0) interface. The superconductor bands are 
modeled by equations (|10[) - (|11[ ). The normal conductance is 
shown for comparison. 
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FIG. 7: The conductance g s is independent of the barrier 
strength, Z, at the energy of the ABS. The latter can be 
checked from Fig. [4] This is shown for two transverse mo- 
menta: p y = 0.15-7T (left); p y = 0.05-7T (right). A e = —Ah = 
0.02. 
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FIG. 8: Relative conductance g s /g n for a disordered inter- 
face, at two different values of transverse momentum. = 
-0.01, A e = 0.025. 
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FIG. 9: Behavior of the relative conductance for different 
disorder values, Z. Destructive interference effects cause a 
zero in the conductance. A e = —Ah = 0.02. 



Since the system is two-dimensional we have to inte- 
grate over the possible incident angles, or over p y . The 
integrated relative conductance is shown in Figure fTUl for 
three interface disorder strengths. The peak structure for 
large values of the barrier strength reveals the existence 
of ABS's. 

In a recent preprint^, a theory is provided that quali- 
tatively predicts the interference effects in the multiband 
superconductor, namely the suppression of conductance 
and the appearance of ABS, in agreement with our find- 
ings. In Re^, the wave function in the superconductor 
is written in the same form as equation (|15|) but it is as- 
sumed that the ratios of the amplitudes E/C and F/D 
are equal (to a phenomenological parameter a introduced 
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FIG. 10: Behavior of the integrated relative differential con- 
ductance, gs/un, for different disorder values, Z. A e = 
-A h = 0.02. 



in Ref 2 ^). Using our waveguide theory for the interface 
matching conditions, we find that the ratios E/C and 
F/D are different, as can be seen from figure fTTI 
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FIG. 11: The ratio of the moduli of amplitudes \E/C\ and 
| F/D | as a function of energy. The transverse momentum 
p v — 0.057T (left panel) and p y = 0.157T (right panel). A e = 
-A h = 0.02. 



V. SUMMARY 

We have introduced a generalization of the quantum 
waveguide theory to determine the appropriate boundary 
conditions for the wave function at the interface between 
a normal metal and a multiband superconductor. We 
have shown that resonant transmission and destructive 
interference effects occur in the sign-reversed scenario for 
pnictide superconductors. Unlike other unconventional 
superconductors, Andreev bound states at finite energies 
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are brought about by these interference effects. 

On the experimental side, polycrystalline samples have 
been used so far. The results obtained above describe an 
interface parallel to the nearest Fe-Fe bonding. There- 
fore, experiments with single crystals are highly desir- 
able. If the edge of the sample is such that the conserva- 
tion of the transverse momentum p y intercetps only one 
FS pocket, existing one band theories apply. The above 
quantum waveguide theory can in principle be used to de- 



scribe other MBS, such as the heavy- fermion materials 22 . 
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